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In quintessence models a scalar field couples with the dominant constituent 
and only acts like a cosmological constant after the onset of the matter dom- 
inated epoch. A generic feature of such solutions, however, is the possibility 
of significant energy density in the scalar field during the radiation dominated 
epoch. This possibility is even greater if the quintessence field begins in a 
kinetic-dominated regime, for example as might be generated at the end of 
"quintessential inflation." As such, these models can be constrained by pri- 
mordial nucleosynthesis and the epoch of photon decoupling. Here, we analyze 
both kinetic dominated and power-law quintessence fields (with and without 
a supergravity correction). We quantify the allowed initial conditions and 
effective-potential parameters. We also deduce constraints on the epoch of 
matter creation at the end of quintessential inflation. 
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1 Introduction 



Observations |T], |2] of Type la supernovae and the power spectrum of the cosmic microwave 
background, together with complementary observational constraints ||, all indicate that 
the universe is presently accelerating. In one interpretation, this acceleration is due to the 
influence of a dominant dark energy with a negative pressure. The simplest candidate for 
such dark energy is a cosmological constant for which the equation of state is u = Pj p = 
— 1. A second possibility is derived from the so-called "quintessence" models, in which 
the dark energy is the result of a scalar field Q slowly evolving along an effective potential 
V(Q). The equation of state is negative —1 < wq < 0, but not necessarily constant. 

Introducing a quintessence field helps to reconcile the fine tuning problem and the 
cosmic coincidence problem associated with a simple cosmological constant. [fH-||ll|. Spe- 



cific forms of the quintessence effective potential can be chosen such that the field Q 
evolves according to an attractor-like solution. Thus, for a wide variety of initial condi- 
tions, the solutions for Q and Q rapidly approach a common evolutionary track. These 
solutions lead naturally to a cross over from an earlier radiation-dominated solution to 
one in which the universe is dominated by dark energy at late times. Another interesting 
possible feature is that such models might naturally arise during matter creation at the 
end of an earlier "quintessential" inflationary epoch [|T2|]. In this case, the Q field emerges 
in a kinetic-dominated regime at energy densities well above the tracker solution. 

It is not yet clear, however, that these models have altogether solved the fine-tuning 
and cosmic-coincidence problems [14| , and there may be some difficulty in embedding 
quintessence models in string theory Nevertheless, several such tracker fields have 
been proposed |I(J whose effective potentials may be suggested by particle physics models 
with dynamical symmetry breaking, by nonperturbative effects ||, by generic kinetic 
terms "fc-essence" in an effective action describing the moduli and massless degrees of 
freedom in string and supergravity theories ||-|f§, or by static and moving branes in a 
dilaton gravity background JT^]. 

A general feature of all such solutions, however, is the possibility for a significant con- 
tribution of the Q field to the total energy density during the radiation-dominated epoch 
as well as the present matter-dominated epoch. The yields of primordial nucleosynthesis 
and the power spectrum of the CMB can be strongly affected by this background energy 
density. Therefore, we utilize primordial nucleosynthesis and the CMB power spectrum 
to constrain viable quintessence models. 



2 Quintessence Field 



A variety of quintessence effective potentials JTUJ or fc-essence effective actions P|-[§] can be 
found in the literature. Observational constraints on such quintessence models have been 
of considerable recent interest [O, [17], [Tj| . In this paper, we describe the work presented 



in |]TTJ] on Q-field and/or kinetic-dominated quintessence models. We concentrate on an 
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inverse power law for the Q field as originally proposed by Ratra and Peebles |T9|, V(Q) = 

m (4+ q )q-^ wher6) M and 

a are parameters. The parameter M in these potentials 
is fixed by the condition that VLq = 0.7 at present. Therefore, Pq(0) = 0.7p c (0) = 
5.7/i 2 x 10" 47 GeV 4 , and M « (p Q (0)Q a ) 1/(Q+4) . If Q is presently near the Planck mass 
and a is not too small (say a^4), this implies a reasonable value || for M which resolves 
the fine tuning problem. 

We also consider a modified form of V(Q) as proposed by [20] based upon the con- 
dition that the quintessence fields be part of supergravity models. The rewriting of the 
effective potential in supergravity depends upon choices of the Kahler potential Jl4|]. The 
flat Kahler potential yields an extra factor of exp{3Q 2 /2m 2 l } ||20| . This comes about by 
imposing the condition that the expectation value of the superpotential vanishes. The 
Ratra potential thus becomes Vsugra{Q) ~~ > M^ A+a ^ Q~ a x exp(3Q 2 /2m 2 l ), where the ex- 
ponential correction becomes largest near the present time as Q — > m p i. This supergravity 
motivated effective potential is known as the SUGRA potential. The fact that this po- 



tential has a minimum for Q = ya/3m p i changes the dynamics. It causes the present 
value of Wq to evolve to a cosmo logical constant (wq ~ —1) much quicker than for the 
bare power-law potential |10| . 

The quintessence field Q obeys the equation of motion Q + 3HQ + dV/dQ = 0, 
where the Hubble parameter H is given from the solution to the Friedmann equation, 
H 2 = (a/a) 2 = l/m 2 pl (p B + p Q ), where m v i = (SnG/3)- 1 / 2 = 4.2 x 10 18 GeV, p B is the 
energy density in background radiation and matter, and a is the cosmic scale factor. As 
the Q field evolves, its contribution to the energy density is given by pq = Q 2 /2 + V(Q). 
Similarly, the pressure in the Q field is Pq = Q 2 /2 — V(Q). The equation of state 
parameter wq is a time-dependent quantity, wq = Pq/pq = 1 — 2V(Q)/ pq, where the 
time dependence derives from the evolution of V(Q) and p(Q). 



3 Nucleosynthesis Constraint 

The quintessence initial conditions are probably set in place near the inflation epoch. 
By the time of the big bang nucleosynthesis (BBN) epoch, many of the possible initial 
conditions will have already achieved the tracker solution. However, for initial conditions 
sufficiently removed from the tracker solution, it is possible that the tracker solution has 
not yet been achieved by the time of BBN at 0.01<T<1 MeV, 10 8 <z <10 10 . 

For many possible initial conditions the tracker solution is obtained before nucleosyn- 
thesis. Along the tracker solution, pq diminishes in a slightly different way than the 
radiation-dominant background energy density. For example, as long as pq « p B , the 
Q-field decays as pq oc a~ 3 ^ 1+WQ \ with wq = (aws — 2)/(a + 2) < %. The equation 
of state wq is only equal to the background equation of state u>b in the limit a — > oo. 
Nevertheless, the tracker solution does not deviate much from pg, even at high redshift for 
most values of a. Hence, one can characterize the nucleosynthesis results by the (nearly 
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constant) ratio Pq/pb during the BBN epoch. If the energy density in the tracker solu- 
tion is close to the background energy density, the nucleosynthesis will be affected by the 
increased expansion rate from the increased total energy density. Such a situation occurs 
for large values of the power-law exponent a. 

A second possibility is that the energy density pq could exceed the tracker solution 
and be comparable to or greater than the background energy density during primordial 
nucleosynthesis. The kinetic energy in the Q field then dominates over the potential 
energy contribution to pq and wq = +1 so that the kinetic energy density diminishes as 
a~ 6 . In this case there could be a significant contribution from pQ during nucleosynthesis 
as the Q field approaches the tracker solution. The strongest constraints on this case 
would arise when pq is comparable to the background energy density near the time of the 
weak-reaction freese out, while the later nuclear- reaction epoch might be unaffected. This 
case is particularly interesting as this kinetic-dominated evolution could be generated by 
an earlier quintessential inflation epoch |12| as described below. 

A final possibility might occur if the Q field approaches the tracker solution from 
below. In this case, the tracker solution may be achieved after the the BBN epoch so that 
a small pQ during BBN is easily guaranteed. However, the ultimate tracker curve might 
still have a large energy density at the later CMB epoch as described below. 

Adding energy density from the Q field tends to increase the universal expansion 
rate. Consequently, the weak reaction rates freeze out at a higher temperature T w . This 
fixes the neutron to proton ratio (n/p fa exp[(m p — m n )/T w \) at a larger value. Since 
most of the free neutrons are converted into He 4 , the primordial helium production is 
increased. Also, since the epoch of nuclear reactions is shortened, the efficiency of burning 
deuterium into helium is diminished and the deuterium abundance also increases. Hence, 
very little excess energy density from the Q field is allowed. The primordial light-element 
abundances deduced from observations have been reviewed by a number of recent papers 



2I|-[23]. There are several outstanding uncertainties. For our purposes |Tl| we adopt 
0.226 < Y p < 0.247 and 2.9 x 10" 5 < D/H < 4.0 x 10" 5 as the most relevant constraints. 

Figure [T] summarizes allowed values of Pq/ Pb at T =1 MeV (z ~ 10 10 ) based upon the 
nucleosynthesis constraints. The region to the right of the curve labeled D/H corresponds 
to models in which the primordial deuterium constraint is satisfied, D/H < 4.0 x 10 -5 . 
The region below the line labeled Y p corresponds to Y p < 0.247. The hatched region 
summarizes allowed values for the energy density in the quintessence field during the 
nucleosynthesis epoch. 

In the present work we deduce an absolute upper limit of 5.6% of the background 
radiation energy density allowed in the quintessence field. This maximum contribution 
is only allowed for 7710 ~ 4.75 or Vt^h 2 « 0.017. A smaller contribution is allowed for 
other values of r/10. Indeed, this optimum value is 4<r less than the value implied by 
the cosmic deuterium abundance @, || VL b h 2 = 0.020 ± 0.001 (la) (r] 10 = 5.46 ± 0.27). 
The independent determinations of VL b h 2 from high-resolution measurements of the power 
spectrum of fluctuations in the cosmic microwave background favor a value even higher. 
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Figure 1: Constraints on the ratio of the energy density in the quintessence field to the 
background energy density Pq/pb (at T = 1 MeV) from the primordial abundances as 
indicated. The allowed region corresponds to Y p < 0.247 and D/H < 4.0 x 10 -5 . 



Both the BOOMERANG ^§ and DASI |27| data sets imply Q b h 2 = 0.022±g;gg| (Iff) 
(7710 = 6.00toii)- The deuterium and CMB constraints together demand that 7710 > 5.19 
which would limit the allowed contribution from the Q field to < 2% of the background 
energy density. 

The most restrictive CMB constraint on derives from demanding a flat universe 
(fltot = 1-0), and marginalizing the likelihood function over all parameters with assumed 
Gaussian priors [26] based upon measurements of large-scale structure and Type la super- 
novae. This gives Vl^h 2 = 0.023 ± 0.003 (Iff). If one adopts this as a most extreme case, 
then Qt,h 2 > 0.020. This would correspond to r/10 > 5.46. From Figure 1 this would imply 
a much more stringent constraint that only about 0.1% of the background energy density 
could be contributed by the Q field. Of course, this is only a Iff constraint, and the upper 
limit to Y p is not well enough established to rule out a contribution to the energy density 
at the 0.1% level. We adopt the conservative constraint of 5.6%. Nevertheless, it is of in- 
terest to explore how the quintessence parameters allowed by BBN might improve should 
the constraints from BBN ever be so tightly defined. Therefore, we will consider 0.1% as 
a conceivable limit that demonstrates the sensitivity to BBN. Even the most conservative 
5.6% limit adopted here corresponds to only about half of the energy density allowed in 
[24] for 3 neutrino flavors. 
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4 Equation of State Constraint 



The Q field must behave like dark energy during the present matter-dominated epoch, 
i.e. the equation of state should be sufficiently negative, (wq = Pq/pq < 0) by the 
present time. We adopt observational constraints on wq from |§. They adopted a most 
conservative approach based upon progressively less reliable data sets. Using the most 
reliable current low-redshift and microwave background measurements, the 2a limits are 
— 1 < wq < —0.2. Factoring in the constraint from Type la supernovae reduces the range 
to — 1 < wq < —0.4. This derives from a concordance analysis of models consistent with 
each observational constraint at the 2a level or better. A combined maximum likelihood 
analysis suggests a smaller range of —0.8 < wq < —0.6 for quintessence models which 
follow the tracker solution, though wq ~ —1 is still allowed in models with nearly a 
constant dark energy. We invoke these three possible limit ranges for the present value of 
wq. 



5 CMB Constraint 



There are two effects on the epoch of photon decoupling to be considered. In the case 
where the energy density in the quintessence field is negligible during photon decoupling 
|IIJ the only effect of the dark energy is to modify the angular distance-redshift relation 
[ 25 1 . The existence of dark energy along the look-back to the surface of last scattering 



shifts the acoustic peaks in the CMB power spectrum to smaller angular scales and larger 
Z-values. The amplitude of the first acoustic peak in the power spectrum also increases, 
but not as much for quintessence models as for a true cosmological constant A. The basic 
features of the observed power spectrum [|26| can be fit |I(J with either of the quintessence 
potentials considered here. For our purposes, this look-back constraint is already satisfied 
by demanding that Qq = 0.7 at the present time. 

The second effect occurs if the energy density in the Q field is a significant fraction of 
the background energy during the epoch of photon decoupling. Then it can increase the 
expansion rate and push the I values for the acoustic peaks to larger values and increase 
their amplitude |25j| . Such an effect has been considered by a number of authors in various 
contexts [Z^, |2£| . For our purposes, we adopt a constraint JI7J based upon the latest 
CMB sky maps of the Boomerang |26j and DASI |27|] collaborations. The density in the 
Q field can not exceed Qq < 0.39 during the epoch of photon decoupling. This implies a 
maximum of Pq/pb = ^q/(1 — ^q) ^0.64 during photon decoupling. 



6 Quintessential Inflation 

Another possible constraint arises if the kinetic term dominates at an early stage. In this 
case pQ ~ Q 2 /2 and pq decreases with scale factor as a -6 . At very early times this kinetic 
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regime can be produced by so-called "quintessential inflation" |12| |. In this paradigm 
entropy in the matter fields comes from gravitational particle production at the end of 
inflation. The universe is presumed to exit from inflation directly into a kinetic-dominated 
quintessence regime during which the matter is generated. An unavoidable consequence 
of this process, however, is the generation of gravitational waves along with matter and 



the quanta of the quintessence field [12|, |29|, 30, ^TJ at the end of inflation. 



6.1 Energy in Quanta and Gravity Waves 

The energy density in created particles can be deduced using quantum field theory in 
curved space-time ]T2|, |30|, [HJ, p B « (l/128)N s Hf(z + l/z 1 + l) i (g 1 /g eff (z)) 1 / 3 , where 
Hi and z\ are the expansion factor and redshift at the matter thermalization epoch at the 
end of quintessential inflation, respectively. The factor of 128 in this expression comes 
from the explicit integration of the particle creation terms. 

When the gravitons and quanta of the Q field are formed at the end of inflation, 



one expects [0 the energy density in gravity waves to be twice the energy density in 
the Q-field quanta (because there are two graviton polarization states). In this paradigm 
then, wherever we have deduced a constraint on pq, it should be taken as the sum of three 
different contributions. One is the dark energy from the vacuum expectation value (pq) of 
the Q field; a second is from the fluctuating part psq of the Q field; and a third is from the 
energy density pew in re hc gravity waves. Thus, we have pq — > ((pq) + Psq + Pgw)- The 
energy density in gravity waves and quanta scales like radiation after inflation, Pgw+Psq oc 
a -4 , while the quintessence field vacuum expectation value evolves as (pq) oc a -6 during 
the kinetic dominated epoch. This epoch following inflation lasts until the energy in the 
Q field falls below the energy in background radiation and matter, pq < ps- 

Thus, for the kinetic dominated initial conditions, gravity waves could be an important 
contributor to the excess energy density during nucleosynthesis. The relative contribution 
of gravity waves and quintessence quanta compared to the background matter fields is 
just given by the relative number of degrees of freedom. At the end of inflation, the 
relative fraction of energy density in quanta and gravity waves is given by (T^j (psq + 
Pgw)/pb = 3/iV s , where N s is the number of ordinary scalar fields at the end of inflation. 
In the minimal supersymmetric model N s = 104. Propagating this ratio to the time of 
nucleosynthesis requires another factor of (gn/di) 1 ^ 3 where g n = 10.75 counts the number 
of effective relativistic degrees of freedom just before electron-positron annihilation, and 
gi counts the number of degrees of freedom during matter thermalization after the end 
of inflation. In the minimal standard model gi = 106.75, but in supergravity models this 
increases to ~ 10 3 . 

Combining these factors we have (psq + Pgw)/ Pb < 0.014, during nucleosynthesis. 
Hence, in this paradigm, the allowed values of pq/pB consistent with nucleosynthesis 
could be reduced from a maximum of 0.056 to 0.042, further tightening the constraints 
deduced here. 
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6.2 Gravity- Wave Spectrum 

There has been considerable recent interest p{| |3(J in the spectrum of gravity waves 



produced in the quintessential inflation paradigm. One might expect that the COBE 
constraints on the spectrum also lead to constraints on the Q field. However, we conclude 
below that no significant constraint on the initial pg or effective potential is derived 
from the gravity wave spectrum. On the other hand, the BBN and CMB gravity-wave 
constraints can be used to provide useful constraints on the quintessential inflation epoch 
as we now describe. Our argument is as follows: The logarithmic gravity-wave energy 
spectrum observed at the present time can be defined in terms of a differential closure 
fraction, VL GW (v) = (1 / Pc)(dpcw I 'dY&v) , where the paw is the present energy density in 
relic gravitons and p c (0) = 3Hq/8tiG = H^m^ is the critical density. This spectrum has 



been derived in several recent papers p9L |30|. It is characterized by spikes at low and high 



frequency. The most stringent constraint at present derives from the COBE limit on the 
tensor component of the CMB power spectrum at low multipoles. There is also a weak 



constraint from the binary pulsar |30| and an integral constraint from nucleosynthesis as 
mentioned above. 

For our purposes, the only possible new constraint comes from the COBE limits on 
the tensor component of the CMB power spectrum. The soft branch in the gravity-wave 
spectrum lies in the frequency range between the present horizon uq = 1.1 x 10^ 18 Q\^ 2 h Hz 
and the decoupling frequency ^d ec (0) = 1-65 x lCT 16 fi]^ 2 /i Hz, where we adopt SIm — 0-3 
for the present matter closure fraction. The constraint on the spectrum can be written 



where, f2 7 = 2.6 x 10~ 5 h 2 is the present closure fraction in photons. The number of 
relativistic degrees of freedom at decoupling is gd ec = 3.36. As noted previously, g\ is 
the number of relativistic degrees of freedom after matter thermalization. In the minimal 
standard model is g\ = 106.75. The quantity H\ is the expansion rate at the end of 
inflation. It is simply related to the kinetic energy pq after inflation, Pq(zi) = Hfm^. 
The logarithmic term in Eq. (|IJ) involves the ratio of the present values of the frequency 
u r (0) characteristic of the start of radiation domination (p B = pg at z = z r ), to the 
frequency characteristic of matter thermalization at the end of inflation ^i(O). This ratio 
is just z/ r (0)/Vi(0) = (z r + 1/zi + l) 2 . The identity p B = Pq at z = z r then gives, 
v r jv\ = N s /128(Hi/m p i) 2 (gi/g r ) 1 ^ 3 , where for the cases of interest g r = 10.75. 

Collecting these terms, we can then use Eq. (|IJ) to deduce a constraint on the ex- 
pansion factor at the end of inflation Hf < 1.4 x 10 _11 mp;. For kinetic dominated 
models, Pq(z{) at the end of inflation is simply related to the energy density pq at z, 
Pq( z ) = Pq( z i)[( z + l)/!- 2 ! + I)] 6 - Similarly the background matter energy density scales 
as p B (z) = pBiz^/g.ffiz)) 1 /^ + \)/{z x + l)] 4 . 
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Considering the present energy density in photons and neutrinos, we can find a relation 
between E x and z x : p 70 + p M = 1.1 x 10" 125 m^ = {N s / 12%) Hf{g x / g dec ) 1 ^{zi + 1)~ 4 . We 
then deduce that z\ < 8.4 x 10 25 and there is only a lower limit on Pq/pb given by the 
constraint on Hi. At our initial epoch z = 10 12 , we find, pq(z)/pb{z) ^5.6 x 10~ 18 . 
This limit is not particularly useful because pq field must exceed ps at Z\ in order for 
the gravitational particle production paradigm to work. The implication is then that all 
initial conditions in which the kinetic term dominates over the background energy at z\ 
are allowed in the quintessential inflation scenario. Hence, we conclude that the gravity- 
wave spectrum does not presently constrain the initial pQ or V(Q) in the quintessential 
inflation model. 

However, the limits on pq/pb < 560 derived from the BBN constraints discussed below 
can be used to place a lower limit on the expansion rate at the end of inflation in this model. 
This in turn implies a lower limit on the redshift for the end of quintessential inflation. 
Thus, we have 2.2 x 1(T 21 < Hf/m^ < 1.45 x 10 -11 , and 1.0 x 10 21 < z x < 8.4 x 10 25 in 
the minimal supersymmetric model. 

This implies that quintessential inflation must end at an energy scale somewhere be- 
tween about 10 8 and 10 13 GeV, well below the Planck scale. By similar reasoning one can 
apply this argument to the gravity- wave spectrum from normal inflation as given in |3(J . 
We deduce an upper limit to the epoch of matter thermalization of Hi < 3.1 x 10~ 10 m 2 z 
which implies z\ < 7.3 x 10 28 [(7(2i)/3.36] 1 / 12 . In this case there is no lower limit from 
BBN as there is no Q field present after inflation. 

7 Results and Discussion 

The equations of motion were evolved for a variety of initial Q field strengths and power- 
law parameters a. As initial conditions, the quintessence field was assumed to begin with 
equipartition, i.e. Q 2 /2 = V(Q), and wq = 0. This seems a natural and not particularly 
restrictive choice, since wq quickly evolves toward the kinetic (wq = +1) or the tracker 
solution depending upon whether one begins above or below the tracker curve. 

Constraints on a and the initial value for the Q-field energy density Pq{z) / Pb{ z ) & t 
z = 10 12 were deduced numerically. These are summarized in Figure § for both: (a) 
the bare Ratra power-law potential; and (b) its SUGRA corrected form. Our constraints 
can easily be converted to closure fraction by VLq = (pq/p_b)/(1 + Pq/pb)- For purposes 
of illustration, we have arbitrarily specified initial conditions at z = 10 12 , corresponding 
to T ~ 10 12 K, roughly just after the time of the QCD epoch. At any time the energy 
density p re i{ z ) in relativistic particles is just p re i{z) = p 7 o (3.36/ 'g e ff(z)) 1/3 (z + l) 4 , where 
p 7 o = 2.0 x 10~ 51 GeV 4 , is the present energy density in microwave background photons, 
and we take g e ff(z) = 10.75 between z = 10 12 and the beginning of BBN just before 
electron-positron annihilation (z ~ 10 10 ). 

The envelope of models which obtain a tracker solution by the epoch of nucleosynthesis 
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Figure 2: Contours of allowed values for a and initial Pq/ Pb ( & t z = 10 12 ) from var- 
ious constraints as indicated for (a) the bare power-law potential, and (b) the SUGRA 
corrected potential. Models in which the tracker solution is obtained by the BBN epoch 
are indicated by the upper and lower curves. Values of a to the right of the lines labeled 
wq = —0.6, —0.4, —0.2 on (a) are excluded by the requirement that the present equation 
of state be sufficiently negative. The BBN constraint for a maximum energy density in 
the Q field of 0.1% (dotted line) and 5.6% (solid line) are also indicated. For the SUGRA 
potential (b) all tracker solutions to the right of the region labeled wq = — .8 are allowed. 



are indicated by upper and lower curved lines in Figures 0a and 0b. The general features 
of these constraints are as follows. If the initial energy density in the Q field is too large, 
the tracker solution is not reached by the time of BBN. The Q-field energy density can 
then significantly exceed the background energy during nucleosynthesis. This situation 
corresponds to the excluded regions on the top of Figures 0a and 0b. All solutions 
consistent with the primordial nucleosynthesis constraints are also consistent with our 
adopted CMB Q-energy constraint as also shown at the top of Figures 0a and 0b. 

The excluded regions at the top and bottom of figures 0a and 0b can be easily 
understood analytically. For example, the excluded (no A) region at the bottom of 
these figures reflects the fact that if pq is initially below the value presently required 
by fl\ = 0.7 it can not evolve toward a larger value. Hence, Pq/ Pb < 2.8 x 10~ 44 is ruled 
out for h = 0.7. Similarly, the "Excluded by BBN" region comes from requiring that 
Pq(zbbn) I Pb( z bbn) < 0.056. For this constraint we are only considering cases in which 
the Q field is approaching the tracker solution from above during nucleosynthesis. Hence, 
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it is in the kinetic regime in which pq oc z 6 while the background scales as z 4 . Thus, 
we have p Q (z = 10 12 )/p B (z = 10 12 ) > 0.056(10 12 /10 10 ) 6 ^ 4 = 560 is excluded. By similar 
reasoning, the "Excluded by CMB" region is given by Pq(zcmb) / 'pb(zcmb) > 0.64, or 
p Q ( z = io 12 )/ p B (z = 10 12 ) > 0.64(10.75/3.36) 1/3 (10 12 /10 3 ) 6 - 4 = 9.4 x 10 17 . 

For the bare Ratra power-law potential (Fig. 2a) the main constraint is simply the 
requirement that the equation of state be sufficiently negative by the present time. The 
sensitivity of the allowed power-law exponent to the equation of state is indicated by the 
wq = -0.2, -0.4, and -0.6 lines on Figure [|a. In the present Q-dominated epoch, lines of 
constant wq must be evaluated numerically. The slight slope to these curves comes from 
the fact that V(Q)/p(Q) has not yet reached unity, i.e. there is still some small kinetic 
contribution to p(Q) and the amount of kinetic contribution depends upon a. 

For the bare Ratra power-law potential, tracker solutions with a<20 are allowed if 
wq < —0.2. The allowed values for a reduce to < 9 and < 2 if the more stringent 
-0.4 and -0.6 constraints are adopted. However, if a is too small, say a < 2, then the 
potential parameter M becomes a very small fraction of the Planck mass and the fine 
tuning problem is reintroduced. 

For models in which the tracker solution is obtained by the time of BBN, the potential 
parameters are only constrained if the most conservative equation of state limit (wq < 
—0.2) and most stringent nucleosynthesis constraint (pq/pb < 0.1%) are adopted. On 
the other hand, independently of the equation of state constraint, nucleosynthesis limits a 
large family of possible kinetic-dominated solutions even though they provide the correct 
present dark energy. 

For the SUGRA-corrected Q fields (Fig. 2b), the constraint from wq is greatly relaxed. 
In fact, wq is sufficiently negative (wq < —0.6) for all a < 10 4 . The reason is that 
all tracker solutions have wq m —1. This is because wq decays much faster toward - 
1 for the SUGRA potential. Also, the potential has a finite minimum which is equal 
to the present dark-energy density. The Q field quickly evolves to near the potential 
minimum and has negligible kinetic energy by the present time. Any potential which 
becomes flat at late times gives wq ~ — 1 and the dark energy looks like a cosmo logical 
constant. All SUGRA models which achieve the tracker solution also have a small pq 
during primordial nucleosynthesis. Hence, there is no constraint from nucleosynthesis 
except for those kinetic-dominated models in which the Q field is still far above the 
tracker solution during the nucleosynthesis epoch. 

We do note, however, that if a lower limit of wq > —0.8 is adopted for tracker solutions 
from ||, then only a power law with a>30 is allowed. This makes the SUGRA potential 
the preferred candidate for quintessence. The large a implies values of M close to the 
Planck mass, thus avoiding any fine tuning problem. However, this potential will be 
constrainable by BBN if the light-element constraints become sufficiently precise to limit 
Pq at the 0.1% level. 
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8 Conclusions 



We conclude that for both the bare Rata inverse power-law potential and its SUGRA- 
corrected form, the main constraints for models which achieve the tracker solution by 
the nucleosynthesis epoch is from the requirement that the equation of state parameter 
becomes sufficiently negative by the present epoch. The main constraint from nucle- 
osynthesis is for models which are kinetic dominated ar the time of nucleosynthesis. The 
SUGRA-corrected potential is the least constrained and avoids the fine tuning problem for 
M. Therefore, it may be the preferred candidate for the quintessence field, although BBN 
may eventually limit this possibility. We also note that the constraints considered here 
provide useful constraints on the regime of matter creation at the end of quintessential 
inflation. 
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11127220, 12047233) of the Ministry of Education, Science, Sports, and Culture of Japan, 
and also by DoE Nuclear Theory Grant (DE-FG02-95-ER40934 at UND). 

References 

[1] P. M. Garnavich, et al, Astrophys. J., 509, (1998) 74; Riess, et al., Astron. J, 116, 
(1998) 1009. 

[2] S. Perlmutter, et al., Nature 391, (1998) 51; Astrophys. J., 517, (1998) 565. 

[3] L. Wang, R.R. Caldwell, J.R Ostriker and P. J. Steinhardt, Astrophys. J., 530, (2000) 
17. 

[4] C. Wetterich, Nucl. Phys., B302, (1988) 668. 

[5] I. Zlatev, L. Wang, and P. J. Steinhardt, Phys. Rev. Lett., 82, (1999) 896. 

[6] T. Chiba, T. Okabe, and M. Yamaguchi, Phys. Rev. D62, (2000) 023511. 

[7] C. Armendariz-Picon, V. Mukhanov, and P. J. Steinhardt, Phys. Rev Lett., 85, 
(2000) 4438. 

[8] C. Armendariz-Picon, V. Mukhanov, and P. J. Steinhardt, Phys. Rev D63, (2001) 
103510. 

[9] P. J. Steinhardt, L. Wang and I. Zlatev, Phys. Rev. D59, (1999) 123504. 
[10] P. Brax, J. Martin, and A. Riazuelo, Phys. Rev D62, (2000) 103505. 
[11] M. Yahiro, G. J. Mathews, K. Ichiki, T. Kajino, Phys. Rev. D, (2002) in press. 



11 



P.J.E. Peebles and A. Vilenkin, Phys. Rev. D59, (1999) 063505. 

C. Kolda and D. H. Lyth, Phys. Lett., B458, (1999) 197; S. Weinberg, in Proc. Dark 
Matter 2000, Marina Del Rey, Feb, 2000, |astro-ph/0005265l (2000). 



E. J. Copeland, N. J. Nunes, and F. Rosati, Phys. Rev. D62, (2000) 123503. 



T. Banks, JHEP submitted (2000), |hep-th/0007T46| ; T. Banks and W. Fischler, 



JHEP submitted (2000), |hep-th/0102077| ; S. Hellerman, N. Kaloper, and L. Susskind, 



JHEP, 06, 003 (2001); W. Fischler et al., JHEP, 07, 003 (2001) |hep-th/010418Tt C. 
Kolda and W. Lahneman, [hep-th/ 01 053001 . 

B. Chen and F-L. Lin, Phys. Rev. D65, (2002) 044007. 

R. Bean, S. H. Hansen, and A. Melchiorri, Phys. Rev. D64, (2001) 103508. 

X. Chen, R. J. Scherrer, and G. Steigman, Phys. Rev D63, (2001) 123504. 

B. Ratra and P.J.E. Peebles, Phys. Rev. D37, (1988) 3406; P.J.E. Peebles and B. Ra- 
tra, Astrophys. J. Lett., 325, (1988) L17. 

P. Brax and J. Martin, Phys. Lett., B468, (1999) 40; Phys.Rev. D61 (2000) 103502; 
P. Brax, J. Martin, A. Riazuelo, Phys. Rev. D62, (2000) 103505. 

K.A. Olive, G. Steigman, and T.P.Walker, Phys. Rep. 333, (2000) 389. 

K.M. Nollett and S. Buries, Phys. Rev. D61, (2000) 123505; S. Buries, K.M. Nollett, 
& M.S. Turner, Phys. Rev. D63, (2001) 063512. 

D. Tytler, J. M. O'Meara, N. Suzuki, & D. Lubin Phys. Rep., 333, (2000) 409. 

K. Freese, F. Adams, J. Frieman, and E. Mottola, in Origin and Distribution of the 
Elements, G. J. Mathews, ed. (World Scientific: Singapore) (1988), pp. 97-115. 

W. Hu, D. Scott, N. Sugiyama, and M. White, Phys. Rev. D52, (1995) 5498. 

C. B. Netterfield, et al. (Boomerang Collaboration), Submitted to Astrophys. J. 
(2001), ^stro-ph/0l044^q . 



C. Pryke, et al. (DASI Collaboration), Submitted to Astrophys. J. (2001), |astro- 
ph/010449U . 



S. Hannestad, Phys. Rev. Lett., 85, (2000) 4203. 

A. Riazuelo, and J.-P. Uzan, Phys. Rev. D62, (2000) 083506. 

M. Giovannini, Class. Quant. Grav 16, (1999) 2905; Phys. Rev. D60, (1999) 123511. 
L. H. Ford, Phys. Rev. D35, (1987) 2955. 



12 



